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Hi  AMIIMIC  T 


During  the  months  of  September  and  October  1971,  the  Naval  Research  Laboratoiy  per  ormed 
optical  measurements  at  the  Naval  Undersea  Research  and  Development  Center  tower  off  Mission  Bay  near 
San  Diego,  California.  The  opticai  measurements  consisted  of  oblique  photographs  and  wide-angie  photog^^^ 
raphy  from  which  the  ocean  spectrum  was  analyzed  as  a  function  of  time  using  the  Sea  PJ"lt’^n“'ysi^  ’ 
technique.  In  conjunction  with  measurements  being  made  simultaneously  by  other  investigators,  an  attempt 
was  made  to  experimentally  determine  the  interaction  of  internal  waves  and  wind-driven  surface  waves_J 
report  discusses  the  results  of  the  opticai  analysis  of  the  sea  photographs  taken  at  15-s  intervals.  The  results 
show  temporal  variations  of  the  spectral  amplitude  in  the  range  of  wavelengths  Dorn  about  3  cm  to  1  m . 
These  temporal  variations  exhibit  changes  as  large  as  15  dB  on  occasion.  Because  or  the  superposition  or 
many  effects,  such  as  tides,  magnitude  and  direction  variations  of  the  wind,  internal-wave  amplitude  con¬ 
tamination;  no  unique  signature  was  deduced  Tor  the  interaction  or  internal  and  surface  waves.  There  is 
evidence  for  a  mild  frequency  dependence  of  the  interaction. 
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abstract 


Dunng  the  months  of  ^Imr  ^Ocjob*  «"•  "  "h 

laboratory  performed  optical  measu  .  ^  near  gan  Diego,  California, 

and  Development  Center ^  tower  ofl ^  widc.angle 
The  optical  measurements  consiste  9  analyzed  as  a  function  of 

photography  from  which  the  o<*an  *  technique  In  conjunction  with 
time  using  the  Sea  Photo  »  “‘tempt 

measurements  being  made  ®'mU  i  “  ^interaction  of  internal  waves  and 
was  made  to  expenmentelly dtacuwe.  the  results  of  the  optical 
wind-driven  surface  waves.  Th  P°  ,  intervais  The  results  show 
analysis  of  th' In  the  range  of  wavelengths 
temporal  venation.  '),  lhe exhibit  changes  a.  large 
from  about  3  cm  to  1  m.  these  re  p  of  many  effects,  such 

as  15  dB  on  occasion.  Because  of  th  **  F?he  internal-wave  ampli- 

as  tides,  magnitude  and  direction  was  deduced  for  the  interaction 

is  cvidence  for  a  miid  frequency  depen‘ 

Hpnre  of  the  interaction. 
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OPTICAL  ANALYSIS  FOR  THE  1971  ARPA  TOWER  EXPERIMENT 


INTRODUCTION 

At  moderate  wind  velocities,  the  sea  surface  appears  to  the  eye  as  a  randomly  occurring 
pattern  of  surface  perturbations  of  never  ending  variability.  Therefore,  it  is  evident  that 
statistical  measures  are  necessary  to  describe  it.  Since  a  point-function  description  would 
have  a  probability  of  occurrence  of  zero,  it  would  be  analytically  monsterous.  Fortunately, 
the  conceptually  simple,  Fourier  decomposition  of  the  two-dimensional  surface  results  in 
fairly  simple  spectral  trends  and  is  very  informative.  With  the  advent  of  Sea  Photograph 
Analysis  (SPA),  the  fundamental  functional,  the  two-dimensional  slope,  or  the  tilt-angle 
spectra  has  become  an  observable.  Recently,  a  need  for  such  a  representation  has  arisen 
to  assist  ARPA  in  the  analysis  of  the  coupling  of  surface-wave  phenomena  with  the  passage 
of  internal  waves.  In  the  past  year,  the  Naval  Research  Laboratory  (Electromagnetic  Scat¬ 
tering  Branch  (Code  5270),  Electronic  Division)  has  assisted  in  the  overall  ARPA  program 
by  determining  in  situ  the  behavior  in  the  presence  of  internal  waves  of  that  portion  of  the 
ocean-wave^  spectrum  lying  between  1  m  and  3  cm. 

The  method  used,  sea  photograph  analysis,  uses  oblique  photographs  of  the  sea  sur¬ 
face  to  establish  the  directional  spectrum  in  the  “high”-frequency  portion  of  the  ocean 
energy  spectrum.  This  is  possible,  since  under  reasonably  generalized  conditions  a  single 
photograph  of  the  ocean  surface  contains  sufficient  information  from  which  to  deduce 
the  directional  distribution  of  constituent  sinusoids  which  make  up  the  highly  complex 
ocean  surface.  Such  statistical  information  is  of  prime  importance  in  any  evaluation  of 
the  mechanisms  of  wave-wave  interactions  or  electromagnetic  scatter  from  the  sea.  This 
report  will  update  the  theoretical  bases  and  experimental  techniques  which  combine  to 
provide  remote  access  to  many  of  the  statistical  properties  of  the  ocean-surface  perturba¬ 
tions,  including  the  time-dependent  spectral  variations  which  are  necessary  in  studying 
internal-wave /surface-wave  interaction. 


COORDINATE  SYSTEMS 

To  study  the  features  of  the  ocean  surface  of  importance  to  optical  analysis,  it  is 
first  necessary  to  assign  a  system  of  coordinates.  Figures  1  and  2  illustrate  the  coordinate 
system  chosen  for  this  analysis.  The  rationale  in  choosing  the  partition  into  a  local  coor¬ 
dinate  system  and  one  locating  the  localized  point  in  relation  to  the  sensor  is  largely  the 
simplicity  involved.  In  general  the  technique  will  use  the  local  system  to  determine  the 
parameters  dealing  with  the  properties  of  a  point  on  the  surface  of  the  water  and  the 
larger  system  to  describe  the  variation  expected  throughout  the  field  of  view  of  the  camera. 

Figure  1  illustrates  the  local  coordinate  system.  Here  C  repiesents  the  unit  vector 
from  the  point  under  consideration  toward  the  camera,  intersecting  the  optical  axis  of 
the  camera  at  the  center  of  the  lens,  N  represents  the  normal  to  the  water  at  the  surface 
point,  and  S  is  the  vector  toward  the  sky  in  the  direction  of  specular  reflection.  The 
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local  cartesian  system  referring  to  the  same  surface  point  is  labeled  XYZ.  The  angle  0  is 
between  X  and  the  projection  of  N  onto  the  XY  plane.  For  a  sinusoidal  surface  disturbance, 
\p  is  constant.  The  angle  of  incidence  of  the  sky  ray  with  respect  to  the  surface  normal  is 
8.  The  depression  angle  0  of  the  surface  point  from  the  camera  is  measured  relative  to  the 
horizontal  plane;  $  is  the  complement  of  0.  The  angle  0  is  the  tilt  angle  of  the  surface 
normal  from  the  zenith  and  essentially  is  related  to  the  slope  of  the  surface  at  a  point. 

The  angle  r?  is  the  normal  spherical  angle  relating  the  zenith  to  the  sky  ray.  Finally,  ?  is 
the  angle  between  the  X  axis  and  the  projection  of  the  sky  ray  onto  the  XY  plane. 

Figure  2  illustrates  the  camera  coordinate  system  (the  angles  6  and  0  are  the  basic 
quantities  used  in  describing  points  on  the  film).  Here  6  is  the  angle  between  the  ray  to 
the  surface  point  and  the  ray  to  the  nearest  point  falling  on  the  projection  of  the  optic 
axis  onto  the  XY  plane.  Also  A  is  the  angle  obtained  by  projecting  6  onto  the  XY  plane 
and  is  therefore  the  angle  connected  with  the  cardinal  directions  (north,  east,  south,  and 

west). 


WAVE  VISIBILITY  IN  TWO  DIMENSIONS 

One  of  the  most  important  aspects  of  SPA  is  wave  visibility.  From  the  onset,  simple 
observations  convince  one  that  surface  waves  are  most  visible  at  oblique  incidence.  Further, 
the  usefulness  of  a  horizontally  polarizing  filter  is  also  apparent.  Beyond  this  it  is  a  dif¬ 
ficult  mathematical  exercise  to  obtain  a  unique  and  quantitative  measurement  of  any  but 
the  most  simple  wave  patterns.  The  full  analysis  of  the  intensity  of  light  arriving  at  an 
observer  at  some  height  above  a  water  surface  is  almost  intractable  for  the  two  spatial 
dimensions  and  a  dimension  of  light  intensity.  Fortunately,  many  of  the  salient  features 
can  be  understood  from  the  simpler  case  which  excludes  the  lateral  spatial  dimension. 

The  extension  to  the  full  situation  should  reveal  little  more  science  and  with  the  present 
choice  of  coordinate  systems  presents  no  real  experimental  difficulty.  For  the  purposes 
of  this  theoretical  discussion,  only  the  two-dimensional  case  will  be  considered.  In  Fig.  2 
the  light  reaching  a  camera  at  C  from  the  ocean  surface  consists  of  three  components: 
that  reflected  from  the  sky,  that  upwel'  -g  component  which  is  scattered  from  beneath 
the  surface,  and  that  contributed  by  the  media  between  the  camera  and  the  surface  point 
under  consideration .  Functionally , 

u  «  TL  +  TC  +  s  (!) 

where  u  is  the  light  energy  density,  T  is  the  reflection  coefficient,  L  is  the  sky  luminance, 

T  is  the  transmission  coefficient,  2  is  the  upwelling  luminance,  and  s  is  the  media  scattering 
function.  More  explicitly, 

lu ($,0)  =  T($  +  0)  L($  +  20)  +  T(£  +  0)  2(F(£  +  0)  +  0)  +  s($)  (2) 

L 

in  the  XZ  plane,  where  F  is  the  refraction  function. 

An  ocean  surface  for  whirh  a  spectrum  can  be  accurately  defined  consists  of  a  smooth 
variation  of  the  normal  angle  0  for  adjacent  points  at  the  air-sea  boundary.  Mathematically 
speaking  a  spectrum  is  undefined  at  wind  speeds  beyond  the  onset  of  whitecapping  because 
the  surface  normal  angle  becomes  multivalued.  Realistically,  the  utility  of  an  ocean-surface 
spectrum  diminishes  with  increasing  whitecap  coverage  for  the  high-frequency  range,  but 
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the  white  water  can  fail  to  seriously  affect  spectral  estimates  in  the  range  of  wavelengths 
containing  most  of  the  wave  energy  until  a  very  sigi  ificant  fraction  of  the  surface  is 
obscured  by  the  v/hitecaps.  With  the  present  state 
to  get  good  photographs  in  such  cases  anyway. 


the  art,  it  is  probably  impossible 


An  examination  of  Eq.  (2)  and  its  partial  derivatives  shows  that  highly  discontinuous 
luminance  functions  yield  maximum  contrast  (0  variation  of  u)  at  the  cost  of  overempha- 
sizmg  the  dependence  of  u  on  particular  values  of  0  and  underemphasizing  the  existence 
of  other  slopes.  This  is  observed  with  sun  glitter  and  is  of  use  in  establishing  slope  distri¬ 
butions  but  is  of  little  value  in  ascertainmg  waveheight  distributions.  This  restriction  to 
smooth  and  monotonic  luminance  function  is  partially  removed  if  it  is  possible  to  ignore 

the  very  high  portion  of  the  wave-number  spectrum.  (See  the  section  on  subresolution 
waves.) 


The  light  scattering  s  from  the  intervening  medium  represents  a  major  limitation  to 
sea  photo  analysis.  It  degrades  the  transmission  of  the  information  inalterably  and  in 
the  case  of  clouds  may  limit  the  altitude  from  which  it  is  possible  to  see  the  ocean.  Such 
scattering  is  generally  a  function  of  depression  angle  because  of  the  different  distances 
light  must  travel  to  reach  the  camera  from  points  on  the  surface.  Unfortunately  s  h  not 
a  uruform  function  of  angle  because  of  atmospheric  turbulence  and  nonhomogeneous 
particulate  distributions;  further,  it  may  be  a  rapid  function  of  time.  This  accounts  for 

fm  I™6"/”611111  ^  thai  for  high  angu,ar  resolution  in  camera  systems,  the  resulting 

imagery  always  seems  degraded  over  that  which  would  be  calculated  from  the  system 
characteristics  or  even  lower  altitude  photos  of  the  same  wave  system. 

In  many  cases  of  interest,  however,  it  is  possible  to  expand  the  luminance  functions 
St^fequltions5  ^  Tet^n  °nly  ,ow'order  terms.  With  this  assumption,  the 


cuo-r  o^o  +  Vo  +  SQ> 

cui  =  riLo  +  r  0l,  +  v0  +  TqCj  +  Sj, 


(3a) 

(3b) 


and 


q  ■  n 


CUn  2Z  (PqLn-q  +  TqCn-q)  +  S„,  (3n) 

q  *  0 

where  all  the  functions  u,  T,  L,  T,  E,  and  s  are  considered  to  have  the  form 


u  =  u0  +  u1(t-^0)  +  -  ••  +  un($-  *0)n, 


(4) 


make  it  possible  to  determine  L  and  9. 
reflection  and  transmission  coefficients 
are  minimal. 


The  measurement  of  u  and  knowledge  of  the 
is  sufficient,  provided  the  contributions  from  s 
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It  is  now  possible  to  compute  the  "sensitivity,”  i.e.,  wave  contrast  that  will  exist  on 
a  photograph  of  a  wave  system.  The  derivative  u  is  taken  to  express  the  sensitivity  to 
variations  of  wave  normal  angle.  Then 

CT1  u0  =  P'L  +  2PL'  +  T't  +  TC'(1  +  F'),  (5) 

where  the  prime  represents  derivatives  with  respect  to  the  argument.  In  addition, 

C“»u{  =  T'L  +  TL'  +  T'6  +  TK'F'  +  s{; 


then 


C  1  u^  =  C_1  u{  +  TL'  +  TC'  -  s(. 


(6) 


(7) 


Thus  knowledge  of  u  “  u(f)  and  of  expressions  for  L'  and  C'  allows  the  evaluation  of  the 
sensitivity  to  wave  variations,  again  with  s  =  0.  Perhaps  a  more  useful  form  of  Eq.  (7)  is 


d  „  v  d  ,  TL'  +  TC'-s. 

—  (In  u)  =  —  (In  u  +  - - i. 

d0  d{  C-l  u 


(8) 


SUBRESOLUTION  WAVES 


Waves  on  the  ocean  surface  at  a  scale  smaller  than  the  resolution  allowed  by  camera- 
film-air  turbidity  effects  can  introduce  errors  in  the  determination  of  L  and  t.  To  illus- 
trate  this,  one  may  consider  Fig.  3.  The  triangular  structure,  when  viewed  at  this  or  more 
shallow  angles,  causes  shadowing.  Even  at  steeper  angles  there  is  still  a  shortening  of  the 
nearly  shadowed  leg.  Obviously,  water-wave  slopes  which  present  a  broadside  aspect 
weigh  more  heavily  on  the  surface  brightness  than  do  slopes  tilted  away  from  the  observer 
and  partially  shadowed  (projecting  smaller  areas  in  the  camera  direction).  Thus  three 
e  ects  arise.  The  first  causes  a  distortion  of  any  sinusoid,  even  resolved  ones.  The  second 
!"^ce°  a  bias  m  the  apparent  distribution  of  normal  angles;  this  bias  is  a  function  of 
both  the  depression  angle  and  the  unresolved  slope  distribution.  (This  will  lead  to  an 
error  in  the  computation  of  L  and  6  and  thereby  affect  the  measured  spectrum.)  The 
third  effect  is  beneficial  in  that  it  may  extend  the  utility  of  SPA.  On  days  for  which  the 
sky  does  not  exhibit  a  "domed”  variation  of  light  intensity,  e.g.,  have  partially  broken 
clouds,  the  subresolution  waves  will  tend  to  smooth  the  apparent  sky  variations.  Small 
waves  tend  to  cover  the  surface  uniformly,  and  consequently  the  apparent  brightness  of 
some  point  on  a  long  wave  is  smeared  over  the  sky.  Although  this  may  tend  to  wash  out 
the  wave  visibility,  it  can  also  tend  to  remove  discontinuities  in  the  sky  brightness. 


Fig,  3  —  Perspective  distortion  resulting  from 
the  oblique  incidence  angle  of  the  camera 
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It  is  necessary  to  use  an  iterative  procedure  to  determine  L  and  E.  This  involves 
obtaining  a  first  estimate  of  L  and  E  and  then  computing  the  spectrum.  The  effective 
surface  tilt  for  the  subresolution  waves  is  then  estimated,  and  new  values  of  L  and  E  are 
obtained.  The  process  is  continued  unl'.l  the  induced  tilt  is  not  affected  by  the  iteration. 
This  effect  should  result  in  second-order  errors  in  the  wave  sensitivity  equation  and  can 
usually  be  ignored. 


OPTICAL  ANALYSIS 

The  optical  analysis  uses  the  relationship  between  light  amplitudes  in  the  two  focal 
planes  of  a  lens.  A  lens  has  a  transfer  function  between  these  planes  which  becomes  a 
pure  Fourier  transform  when  used  with  collimated,  monochromatic  light.  When  a  positive 
‘  ransparency  of  a  scene  is  inserted  into  the  collimated  beam,  an  image  is  reconstructed 
(just  after  tlie  transparency)  which  to  within  the  camera-film  limitations  is  exactly  that 
which  existed  at  the  camera  due  to  the  scone.  In  water-wave  photography,  the  small 
amplitudes  of  the  waves  (an  inherent  characteristic  of  water  waves  in  the  gravity  range) 
provide  a  linear  relation  between  amplitudes  and  the  surface  normal  angle.  Therefore  the 
light  amplitude  distribution  in  the  other  focal  plane  is  the  Fourier  transform  of  the  surface- 
normal-angle  variations,  and  the  light  intensity  is  proportional  to  the  energy  spectrum  of 
the  water  waves. 
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Fig.  4  —  Optical  analysii  system 


Figure  4  is  a  sketch  of  the  usual  optical  system.  This  system  does  not  use  collimated 
laser  beam  because  of  the  extra  lenses  required  and  the  little  benefit  it  serves.  Instead  the 
lens  images  the  spatial-filter  pinhole  at  the  Fourier  plane;  the  effective  focal  length  for  the 
transform  is  just  the  distance  from  the  transparency  to  the  Fourier  plane.  This  allows  the 
arraignment  to  have  more  versatility  in  that  differing  illumination  areas  and  effective  focal 
lengths  can  be  generated  simply  by  moving  the  transparency  relative  to  the  lens. 


PHOTOGRAPHIC  SPECTRUM 
Goodman  (1)  writes  that 

2 

(9) 

in  which  t(x,  y)  is  the  amplitude  transmittance  function  of  the  scene  for  a  uniform  distri¬ 
bution  of  amplitude  A,  X  is  the  optical  wavelength,  f  is  the  focal  length,  x  and  y  are  the 
distances  in  the  scene  plane,  and  xf  and  y^  are  the  distances  in  the  Fourier  plane.  The 
integral  is  equivalent  to  a  Fourier  transform. 


If(Xf'  yt)  =  M2  [_  J_  t(x* y)  exp  [“*  |?{xxf +  yyr}]  dxdy 
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For  SPA  the  transmittance  function  associated  with  the  film  is 

t(x,  y)=  [l  +  0(x,  y)J  t0.  (10) 

The  laser  beam  is  not  a  uniform  distribution  but  has  a  Gaussian  intensity  distribution,  so 
the  weighted  amplitude  transmitted  becomes 

At(x,  y)  =  A0to  Jj  y)]  exp  ~  ] '  (11) 


Consequently 


i,-(A?)0)1  3{e-o>!  *  y2i/»’} 

(12) 

where  0  denotes  a  Fourier  transform.  The  second  term  in  Eq.  (12)  reduces  to  a  convolution 
of  the  spectra,  which  produces  a  smoothing  of  the  0  spectrum.  Using  the  Four;<;r-transform 
identities 

F(s)  = 

J  f(x)  e  l2nxt  dx, 

—  m 

(13a) 

f(x)  - 

j  F(s)  ei2"**  ds, 

(13b) 

3{f(ax)}  -  jij-  F(|), 

(13c) 

and 

gjg-irx2}  =  e'*’'*2, 

(13d) 

we  arrive  at  an  explicit  form  of  Eq.  (12): 

-(*£) 

1  [lGl^^0,0.*)t(i^)’ 

(14) 

where 

3  {e"(x*  + 

y*)/o2}=  n02  eXp[  (X2+y2)j  =  G, 

(15a) 

^{0(x,  y)}=  0, 

(15b) 

and  the  symbol  *  denotes  a  convolution. 


The  middle  term  of  Eq.  (14)  is  weighted  so  heavily  by  the  Gaussian  term  that  it  is 
to  be  ignored,  as  one  does  the  first  term,  outside  a  region  surrounding  the  origin.  The 
information  is  suppressed  near  the  origin  due  to  the  strongly  peaked  Gaussian  function. 
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This  can  be  looked  on  as  a  physical  manifestation  of  sampling  theory,  i.e.,  the  under- 
sampled  wavelengths  are  obscured  by  the  dc  term.  Then 

,  [ xf  yf\ 

{kx,  ky}  =  27r\\f  *  XfJ’ 

where  G2  is  the  dc  term  and  kx  and  ky  refer  to  film  wavenumbers,  and  letting 


one  obtains 


U  2  t  - 

<t>  »  4ir2o4  —  7 -  =  ^  *  G 

V  i*2  l0 


u2  lo 
0 


The  energy  spectrum  is  related  to  the  square  of  the  Fourier  tmnsform  d2  by 


th+Ak  _ 

1 0(1 

.  4ffzo4  uo  r 

*0  %2  Jk 


1 0  (k)  1 2  dk 


The  intensity  in  the  Fourier-transform  plane  is  constant  over  a  region  inversely  related  to 
rrr^ure  2/2  o  Of  the  Gaussian  illumination  in  the  scene  plane-  Then 


Ak„  =  Akv  = 


-  2tt 


2 ^2  o' 


For  the  circular  region, 


/A k\2  _  if3 

ffT 


and 


Equation  (23)  is  valid  for  the  region  outside  the  origin, 
and  the  photographic  spectrum  is  an  observable. 


(23) 

All  the  quantities  are  measurable, 
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DIRECTIONAL  DISTORTION 


The  snapping  from  real  space  to  film  space  induces  a  distortion,  due  to  the  obliquity. 
Thus  letting  M  represent  a  mapping  operator,  one  has 


where 


M  *(X,Y)  -  t(x,y) 


D 

F 


KX 


esc  0,  and  ky 


(24) 


where  D  is  the  distance  from  the  camera  to  the  surface  point  and  F  is  the  focal  length  of 
the  camera.  As  a  result,  the  photographic  transform  will  exhibit  a  more  restricted  directional 
distribution  than  that  actually  existing  on  the  ocean  surface. 


The  effect  of  oblique  perspective  only  complicates  the  problem  of  obtaining  the  full 
directional  spectra.  To  determine  the  angle  ft  in  the  photo  spectrum  related  to  a  wave  on 
the  ocean  at  a  direction  i p,  consider 


K  =  i  cos  i|/  +  j  sin 

(25a) 

■“('£*♦  I**)* 

(25b) 

k  cos  \p 

.  .  -  cos  ft  ~  > 

Ikl  Vcos2  ip  +  sin20sin20 

(25c) 

and 

ft  =  cos'1  [1  +  tan2  \p  sin2  O]1^2  (26) 

Examination  of  Eq.  (26)  shows  that  obliquity  narrows  the  directional  characteristics  for 
almost  all  \ p . 


OCEAN  SPECTRUM 


It  is  necessary  to  establish  a  relationship  between  the  photographic 
obtained  and  the  real  ocean  spectrum.  The  connection  is  provided  by 

spectrum  thus  far 

3  {«**>}-  B: p(i)- 

(13c) 

where 

3  {f  (x)  }  =  F  (s). 

(27) 

Extending  this  to  two  dimensions  and  operating  on  the  Fourier  transform 

3  {0(X,Y)}=  0 (Kx,Ky) 

(28) 

yields 

3  ^  ky). 

(29) 
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Equations  (28)  and  (29)  represent  the  same  quantity,  the  0  coefficients  being  in  effect  only 
a  change  in  the  dimensions  used  to  measure  X  and  Y.  Therefore, 

^ocean  (Kx,  Ky)  -0x0y  ^pho*,o 

where  the  word  subscripts  are  used  to  keep  track  of  the  “space”  in  which  the  quantity 
is  defined.  Invoking  Eq.  (19)  once  again  yields 

«K  +  AK 

4>0(KX,KY)  =  |0o  (Kx,Ky)|2  dKxdKY  (31) 

JK 


and 

^ocean  (KX,  Ky)  “My  ^photo  (kx*  ky)- 

In  this  case  the  0  *s  are  the  factors  relating  the  scales  on  the  sea  and  photo  and  are  given  by 

{**’M  =  FsinO  {sto?*  *}’  (33) 

hence  one  obtains 

^ocean  (*X.  Ky)  csC3  0  *  photo  (j  esc2  0  Kx,^CSC0KYy  (34) 


LONG-WAVELENGTH  SPECTRUM 

Optical  analysis  performs  exactly  and  very  rapidly  a  mathematical  transformation 
which  is  otherwise  lengthly  and  somewhat  difficult.  The  limitations  of  SPA  basically  reside 
in  the  difficulty  encountered  in  recording  an  appropriate  photograph.  The  long-wavelength 
range  of  the  water-wave  spectra  requires  a  photo  taken  from  an  altitude  sufficiently  high 
such  that  the  transform  falls  outside  the  dc  term.  From  such  altitudes,  many  difficulties 
may  arise  which  would  preclude  good  photography,  i.e.,  haze  and  cloud  coverage.  Wave 
photos  of  wide  angular  coverage  from  lower  altitudes  would  have  too  much  perspective 
distortion  to  yield  reliable  optical  spectra.  It  is  possible,  however,  to  exact  such  low- 
frequency  information  from  a  single  photograph  by  analyzing  the  intensity  variations  along 
a  line  through  the  reconstructed  image  and  by  removing  effects  of  differing  perspective. 

Then  if  a  number  of  such  photos  are  analyzed  and  the  spectra  are  averaged,  the  long-wave¬ 
length  spectral  components  can  be  measured  with  good  statistical  accuracy. 

The  method  used  is  to  scan  the  photograph  with  either  a  microdensitometer  or  a 
flying-spot  scanner,  to  digitize  the  results,  and  to  perform  the  spectral  calculations  by 
computer.  A  new  scan  direction  is  required  for  each  wave  direction.  Either  the  proper 
scan  direction  can  be  obtained  from  the  directivity  of  the  optical  spectra,  or  the  synoptic 
situation  can  be  used  to  infer  the  proper  direction. 

It  is  possible  to  approximate  optically  the  technique  of  averaging  spectra.  The  approach 
makes  a  square  matrix  or  mosaic  of  individual  photos  and  rephotographs  the  composite  to 
obtain  a  single  photo  of  appropriate  size.  The  resulting  spectrum  contains  lower  frequency 
information.  This  technique  creates  about  as  many  problems  as  it  solves  and  to  date  has 
proved  decidedly  inferior  to  the  scanning  technique. 
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Pa,  of  the  great  utility  o,  the  utould^  £S  su" 

SrSSLJSSStSJS  necessary  cognations  is  significant*  " 
scanning  and  digital-computing  techniques  can  again  be  competltiv  . 


SLOPE,  WAVEHE1GHT,  AND  FREQUENCY 
SPECTRA 

The  techniques  tor  obtaining  a  rea.  ~ 

spectra  are,  however,  rather  more  unami  f  iliar  territory  (first  to  slope  spectra,  then 
heights.  Making  the  hransformatio  ctra)  involves  assumptions  which  degrade 

to  waveheight  spectra,  anu  finally  t  q  y  .  .  •  tbat  for  water  waves  the 

the  confidence  in  the  final  results.  The  first  step  is  wa^  .fi 

tangent  of  the  wave  normal  is  a^ost  aiway^  ne  y  q  \  do  occur.  Carrying  the 

spite  of  the  fact  that  even  at  moderate  wind  velocitie  steep  P  difficulties  *hich 

rnSS  Sn"  5MSSS  dSrEuthor  has  lound  ,4,  that 

*(«;K,*)-K**(h;K,*).  (35) 

This  notation,  (<f> ;  K,  ^ ),  is  used  to  denote again  ^"obtain  frequency 
and  azimuth  angle  It  is  now  necessary  ‘  g  in  the  ocean.  This  introduces 

St; rf“^eSTthe  only  approximate  equality  existing  in  the 
classic  equation 


u2  =  gK, 


(36) 


where  ci  is  the  angular  frequency  o,  ^"^""^ce  Twh"  T 
only  for  a  single  wavelength  propagating  on  am n ‘f? Iv  smSand  probably  do  not 

no  currents.  Again,  however,  the  to™*™ ooti^aHy  (Recent  measurements  in  the 
influence  the  first-order  spectra  to  be  obtained  optica  J’ (  nrobablv  due  to  surface  currents 
ultracravity  region  do  exhibit  first-order  spectral  variations,  probat y  due ^to 
S  by  the  orbital  motion  of  .well  waves.)  Withm  these  hm.tat.ons  (2) 


(h;cj) 


2cu_3  f2" 

'  82  Jo 


$  (h;  K,  <l<) 


(37) 


wave  poles  and  can  be  used  to  evaluate  optical  measurements. 

The  distinction  in  the  basic  data  does  nwke  X^U^aTng^S'oTer 

the  different  measurement  systems 

quite  differently  . 
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EMPIRICAL  SPECTRA 

There  exist  empirical  frequency  spectral  forms  obtained  from  many  observations  of 
“fully  developed”  seas  (3).  These  are  one-parameter  distributions  relating  wind  velocity  to 
the  spectrum  which  results  when  equilibrium  is  established.  Just  as  temperature  describes 
thermal  equilibrium,  wind  velocity  seems  to  be  adequate  for  ocean  spectra  at  least  at  long 
wavelengths.  Transient  states  are,  however,  many  parameter  states  and  present  a  significantly 
more  complex  problem. 

A  modern  spectral  form  is  due  to  Pierson-Moskowitz  (3)  and  has  the  form 

<J»(h;co)=^exp{-7(^-)4}  (38) 

a  =  8.10  X  10‘3,  7  =  0.74,  andco0  =  gU^g  B,  where  U19  5  is  the  wind  velocity  at  19.5  m 
above  the  mean  sea  level. 


As  a  consequence  of  this  spectral  form,  many  of  the  oceanographic  observables  are 
specified  in  terms  of  the  wind  velocity.  Among  these  are 


Hi/3  =  2.14  X  10-2  U,9  6, 

(39) 

um«  ”  0.877  U,1,  ,, 

(40) 

*0.81  ^  U„.s, 

(41) 

and  i 

Hrms  =  ^  Hl/3 

(42) 

Hj/3  is  the  average  waveheight  of  the  1/3  largest  waves;  wmax  is  the  frequency  associated 
with  the  peak  in  the  spectrum;  Tave  is  the  average  period  of  the  predominant  wave;  Hrm, 

is  the  rms  waveheight. 

Using  Eqs.  (37)  and  (38),  it  is  possible  to  express  the  directional  spectrum  in  the 

form 


4>  (h ;  K,  iJO  =  K'4  P  (K,  \p ),  (43) 

in  the  region  for  which  the  exponential  factor  is  constant,  where  P(K,  ip)  is  a  dimensionless 
“pattern”  of  the  spectrum,  for  which  the  integral  from  0  to  2tt  equals  the  constant  a  in 
this  region.  Thus  optically  the  function  usually  plotted  is 

K2  <!>(*;  K,  i>)  =  K4  <t>  (h;  K,*)  =  P(K,*).  (44) 

This  removes  the  strong  trend  in  the  spectrum  as  wavenumber  varies.  Further,  it  is  easy 
to  determine  the  onset  of  the  expontential  factor,  enabling  one  to  assign  a  wind  velocity 
to  a  given  spectrum,  assuming  equilibrium  is  established. 
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Substituting  the  spectral  form  of  Eq.  (43)  into  Eq.  (37)  and  using 
spectrum  on  the  left,  one  obtains 


ag^ 


exp 


T(“Ui9.e) 


2cj3 

g2 


K":  P(K) 


using,  to2  =  gK  and  ex  =  1  —  x,  x  <<  1,  yields 

g4 


1-7 


““  Uj9.5 


=  2P(K). 


the  Pierson-Moskowitz 

(45) 

(46) 


Using  full  optical  analysis,  the  right  side  can  be  calculated  to  give  the  value  of  k  for  which 
some  fraction  of  saturation  is  reached,  and  the  left  side  brackets  must  be  the  same  fraction. 
Thus  Ujg  5  is  estimated. 


TRANSIENT  STATES:  STATIONARITY  AND 
HOMOGENEITY 

The  real  wind  over  a  region  of  deep  ocean  varies  with  time  and  space,  whereas  the 
major  dissipation  and  propagation  characteristics  are  largely  invariant  (ignoring  boundaries, 
oil  slicks,  etc.).  Thus  observed  spectra  are  often  a  superposition  of  spectra  generated  at 
different  times  and  places.  Usually  a  given  area  is  designated  as  having  “sea,”  i.e.,  a  portion 
of  the  spectrum  resulting  from  local  winds,  and  “swell,”  i.e.,  spectra  resulting  from  past 
storms  at  other  locales.  The  wave  dissipation  per  unit  length  can  be  thought  of  as  being 
almost  a  constant;  thus  short  wavelengths  damp  out  more  rapidly  than  longer  wavelengths. 

A  spectrum  can  exhibit  a  number  of  peaks,  such  as  an  ultra-gravity-wave  peak  due  to 
inhomogeneous  values  of  the  wind  velocity,  a  gravity -wave  peak  due  to  the  stationary 
component  of  the  local  wind,  a  long-gravity-wave  peak  associated  with  some  distant  storm, 
and  perhaps  a  very-long-wave  peak  due  to  tidal  and  other  motions.  There  may  even  be  a 
capillary-wave  peak  at  light  winds  arising  from  the  local  winds  but  caused  by  the  effects 
of  the  surface  tension  altering  the  spectral  form  for  equilibrium.  Such  a  transient  state  is 
very  common,  and  no  one  parameter  description  can  be  adequate  for  a  complete  ocean 
definition. 

Sea  Photo  Analysis  provides  a  technique  whereby  some  of  these  effects  can  be  studied. 
Photographs  can  be  taken  at  one  place  as  a  function  of  time  or  after  short  time  intervals 
over  seperate  parts  of  the  ocean.  In  this  way  stationarity  and  homogeneity  can  be  evaluated 
and  the  effects  on  other  sensors  estimated. 


TOWER  INSTALLATION 

The  experimental  manifestation  of  the  preceeding  theoretical  considerations  at  the 
Naval  Undersea  Research  and  Development  Center  (NURDC)  tower  in  San  Diego  was  as 
follows.  A  35-mm-format  Nikon  camera  with  a  50-mm  lens,  a  horizontal  polarizing  filter, 
and  a  250-exposure  back  was  detachably  mounted  to  the  tower  at  a  point  62  ft  above 
mean  sea  level.  The  camera  mount  was  affixed  to  the  wench  driven  sled  to  allow  the 
camera  to  be  raised  or  lowered  as  needed  for  film  changes  or  demounting.  The  mount  was 
designed  to  secure  the  camera  in  a  40°  depression  angle  from  the  horizon  with  the  35-mm-film 
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Fig.  5  —  Projections  of  the  camera  field-of-view  onto 
the  ocean  surface 

dimension  subtending  depression  angle  and  the  24-mm  axis  covering  azimuth  angle.  The 
optic  axis  normally  was  oriented  in  a  due  west  direction,  although  provisions  were  made 
to  accomodate  a  due  north  look  angle.  Figure  5  shows  the  geometric  configuration. 

Previous  tests  had  shown  that  the  light  variations  were  sufficiently  large  (as  a  function 
of  time)  to  preclude  the  use  of  a  high-contrast,  low-latitude  film.  Consequently  all  the 
photographs  were  taken  using  Panatomic  X  film,  which  has  an  ASA  rating  of  about  32 
and  a  gamma  of  less  than  1.  Exposures  normally  ran  1/125  with  an  aiperture  of  f/2.8. 
Photographs  of  the  ocean  surface  were  taken  every  15  s.  This  allowed  the  250-exposure 
capacity  to  provide  a  data  run  matching  the  normal  64-min  radar  run. 

Although  it  was  not  mandatory  to  obtain  absolutely  calibrated  spectra,  fish-eye 
photographs  were  taken  at  various  intervals  during  each  day.  This  imagery  was  intended 
to  allow  absolute  calibration  of  the  spectral  data  should  that  be  desired  m  the  future. 

These  photographs  were  taken  by  pointing  the  hand  held  camera  toward  the  horizon  in 
the  direction  the  50  mm  camera  was  aimed.  The  180  or  hemispherical  field  of  view 
allows  the  relative  luminance  of  the  sky  and  the  water  to  be  determined  directly  improving 
the  accuracy  of  the  Eqs.  (7)  and  (8)  used  in  calculating  wave  sensitivity.  This  also  allows 
a  direct  measurement  of  the  subresolution  rms  slope. 
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DATA  ANALYSIS 

The  optical  data  analysis  uses  a  tensjrith  an  ,^„aLl  “P''tute 

The  spot  s,e  in  this  case  ,s 

5.5  mm  in  diameter. 

The  normal  means  of  reading  out  the  t  d^Sti^^S’howewr,  proved  to 

scan  a  photomultiplier  ^^^^the^arge  number  of  photographs  accumulated 
be  far  too  slow  to  conveniently  analyze  the  larg  with  a  significant 

during  the  field  work.  This  is  true  even  when  on* thc*e  d  ^  ^  reduction>  tWo  masks 

overlap  with  the  other  sensors  are  COI\slder  .  ..  ht  attribution  and  sequentially  read 

were  fabricated  which  could  be  moved  through  the  light  interval.  This 

the  integrated  light  intensity  at  some  mean  fixed-percentage  bandwidth  at  a  set 

allowed^  tabulation  of  the  light  intensities  falhng  in  0fthe  masks.  This  mask 

of  discrete  ocean  wavelengths.  Figur ■  hoies  to  transmit  light  at  any 

is  backed  by  a  slit  of  proper  size  to  allow  only  the  ambiguous  other  half.  It 

•—  «  <*«*  ptoportional  to  P(K; w  of  Eq; 

The  use  of  the  masks  was  not  without  *f^y' pj^cTtoleraMK,  however,  provide 
accurately,  especially  for  the  «y  low PJ  small  pinholes  very  new,  the  dc 

teri^show1  e^eme^uctuatTons  is  measured  become^nsignificant 

h  — ths  tM  the  rad” 

used  in  the  radar  system. 


o 


O 


o 


O 


o 


O 


o 


o 


Fig.  6  -  Mask  design  (enlarged  5  times) 
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EXPERIMENTAL  RESULTS 

Figures  7  through  70  represent  the  basic  results  of  the  optical  data  analysis.  Pre¬ 
sented  are  eleven  distinct  runn  on  eight  different  days.  Data  were  taken  on  other  days, 
but  the  analysis  was  sufficiently  lengthy  to  preclude  total  analysis  and  made  it  possible 
to  reduce  only  those  data  having  a  significant  overlap  with  other  sensors.  In  the  suc¬ 
ceeding  paragraphs,  each  data  run  is  discussed  to  point  out  the  significant  features  ot  the 
data  and  to  .ndicate  the  limitations  pertaining  to  the  interpretation  of  the  results.  Most 
of  the  graphs  are  plotted  with  the  negative  sign  of  the  ordinate  suppressed.  Increasing 
spectral  energy  is  therefore  in  the  positive  y  direction.  The  abscissa  is  plotted  as  Green¬ 
wich  Mean  Time  (GMT),  which  differed  by  7  hr  from  the  local  time.  For  all  graphs  the 
spectra  was  sampled  at  15  s  intervals. 

Figures  7  through  10  are  plots  of  the  spectral  variation  for  the  Sept.  16  run  between 
the  times  of  1637  and  1735  GMT.  The  mask  used  averaged  the  spectrum  over  a  ±25% 
band  of  wavelengths.  This  average  operates  over  both  dimensions  in  wavenumber  space. 

It  is  evident  from  these  curves  that  there  are  large  deviations  from  the  mean  of  the  spec¬ 
tral  amplitude  over  this  time  period  with  fine  structured  variations  superimposed.  This  is 
characteristic  of  all  the  runs.  The  large  deviations,  having  an  internal-wave-like  periodicity 
(IW-like),  can  arise  due  to  nonstationary  wind  conditions,  inhomogeneous  wind  field,  or 
internal-wave  (IW)  effects.  The  fine  structure  could  be  due  to  at  least  four  causes,  which 
are  listed  in  the  order  of  importance.  First,  the  folding  of  the  spectra  by  the  existence 
of  swell.  This  effect  is  most  predominant  and  holds  for  all  the  measurements  made  during 
the  experimentation  at  the  tower.  Second,  it  is  possible  that  the  ocean  surface  is  com¬ 
posed  of  numerous  wave  packets  or  patches  of  fairly  monochromatic  energy.  If  the  num¬ 
ber  of  these  packets  in  a  given  analyzed  area  is  small,  the  expected  variation  could  be 
quite  large.  Third,  there  is  some  positioning  error  involved  in  the  mask  placement  from 
frame  to  frame.  As  a  result  there  is  an  imprecision  of  the  data  which  is  most  significant 
for  the  long  wavelengths  but  becomes  minimal  as  the  wavelength  is  shortened.  Finally, 
there  is  a  small  variability  associated  with  the  mean  exposure  level  of  the  photographs. 

This  factor  is  quite  small  and  was  measured  during  the  data  analysis  and  found  to  be 
insignificant.  As  with  all  the  results  of  this  experiment,  the  magnitude  of  thr*  IW-like 
variability  diminishes  toward  the  shorter  wavelengths.  This  is  most  likely  due  to  the 
“frequency”  response,  referred  to  in  optics  as  the  modulation  transfer  function  (MTF), 
of  the  camera-film  system.  The  wave  information  is  of  sufficiently  low  contrast  that 
the  “roll  off”  is  significant  as  the  recorded  film  wavenumber  approaches  20  lines/mm. 

Figures  11  through  14  also  pertain  to  Sept.  16,  but  in  this  case  the  camera  was 
looking  north.  This  orientation,  because  of  wave-visibility  considerations,  predominately 
measures  spectral  components  perpendicular  to  the  internal  waves.  Such  information  is 
of  value  in  analyzing  the  coupling  at  high  angular  differences  between  the  wave  systems. 
Again  there  is  evidence  of  the  fine  structure  associated  primarily  with  the  effects  of  under¬ 
sampled  swell.  There  is,  however,  the  longer  term  fluctuations  probably  due  to  internal 
waves  or  due  to  variations  of  wind  velocity  or  direction.  A  decision  between  the  causes 
would  require  a  cross-correlation  with  the  other  sensors  aboard  the  tower,  such  as  the 
thermocline  fluctuation  and  wind  records.  Compared  with  the  data  taken  during  the 
morning,  the  long-term  variations  are  significantly  smaller,  about  3  dB  compared  to  7  or 
8  dB,  a  result  probably  due  more  to  the  look  angle  than  to  the  time  difference. 
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Figures  15  through  18  are  the  results  of  optically  analyzing  the  photographs  taken 
on  Sept.  28  from  1846  to  1943  GMT.  Here  there  are  two  regions  of  very  strong  suppres¬ 
sion  of  the  wave  spectrum.  At  about  1852  and  1900  GMT,  the  spectral  amplitude  de¬ 
creases  by  a  factor  of  ten.  In  all  probability,  such  extreme  variation  is  the  result  of  a 
“slick”-like  section  of  the  water.  The  width  of  the  slick  is  about  3  to  4  min  in  time, 
probably  corresponding  to  18  to  24  m  of  spatial  extent  if  the  slick  is  moving  at  a  rate 
of  10  cm/s.  Visible  slicks  were  not  normally  this  wide.  It  is  possible  that  the  current- 
wave  interactions  which  tend  to  compact  biological  material  may  have  a  longer  range 
effect  than  that  apparent  to  the  eye.  Again,  correlation  of  these  figures  with  the  slick 
photography  might  prove  to  be  of  value.  The  24-cm  curve  shows  considerable  variability 
of  the  fine  structure.  This  might  result  from  inordinate  positioning  errors  in  the  data 
reduction,  but  this  is  unlikely,  since  the  same  variability  extends  into  the  higher  frequencies 
as  well.  A  more  probable  reason  is  that  a  fairly  significant  swell  existed  on  this  day  which 
caused  severe  changes  in  the  spectrum,  due  to  strong  orbital  currents.  The  fact  that  the 
variability  remains  when  the  analyzed  region  is  many  wavelengths  long  in  both  dimensions 
supports  the  swell  hypothesis. 

Figures  19  through  21  display  the  variation  in  the  spectrum  for  Oct.  1  from  1628 
to  about  1731  GMT.  For  this  run  the  mask  used  averaged  over  only  a  ±5%  region  in 
wavenumber  space.  With  this  more  restricted  frequency-plane  average,  the  fine  structure 
variability  is  increased  significantly.  There  is  a  strong  trend  in  the  data  with  evidence 
for  IW  period  modulation  of  small  magnitude.  The  6-cm  curve  shows  the  small  amplitude 
modulation  clearly,  but  with  such  a  strong  trend,  probably  due  to  a  changing  wind,  this 
fluctuation  may  also  be  due  to  the  temporal  instability  of  the  wind  field. 
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Figures  22  through  29  are  time  histories  of  spectral  fluctuations  for  Oct.  4  from 
1835  to  1928  with  a  ±5%  frequency-plane  average.  This  is  one  of  the  most  unique  sets 
of  data  accumulated  during  the  entire  program.  During  this  time  interval,  the  wind 
record  shows  the  wind  to  be  quite  constant  in  direction  and  magnitude.  During  the  entire 
run,  the  wind  was  270°  t  30°,  and  the  magnitude  was  between  4  and  6  knots,  except  for 
the  interval  1902-1905,  during  which  it  rose  to  about  7  knots.  The  time  histories,  there¬ 
fore,  are  very  closely  those  for  an  equilibrium  condition.  The  graphs  show  that  there 
are  two  slick-like  regions,  one  occurring  at  about  1849  and  another  at  about  1923.  The 
suppression  of  the  ultragravity  waves  at  these  times  is  about  10  dB.  During  the  period 
between  1855  and  1905,  the  wave  spectra  rose  about  3  dB  above  the  equilibrium  level, 
just  before  the  wind  exceeded  the  6-knot  bound  at  1902.  This  region  of  the  internal 
waveform  is  most  likely  a  region  of  differential  roughening  due  to  the  interaction  of  the 
internal -wave/surface-wave  system.  A  further  anomalous  feature  is  the  very  steep  recovery 
after  the  first  slick  at  1850.  Just  to  the  rear  of  the  slick,  the  spectral  amplitude  is  seen  to 
rise  for  a  very  restricted  time  interval  and  then  fall  back  about  3  dB  before  the  overshoot 
starting  at  1856.  The  35-min  period  for  the  arrival  of  internal  waves  is  somewhat  unusual 
but  may  be  accounted  for  by  the  effect  of  tidal  currents  slowing  the  apparent  rate  of 
progression  of  the  wave. 
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Figure  29 


,  .  ,  fv,n  interval  about  the  slick  occurring 

Figures  30  and  31  Me  a  more  ££  aU  the  olher  data.  First,  to  ob- 

at  1849.  Hiese  data  *ere  reduced  quite  drfteenj  ^  superimpolling  the  exposure  ot 
tain  some  smoothing,  the  spectra  ,  e  the  effects  of  swell  or  any  variability  intro- 

four  consecutive  frames.  This  sh°“ld  .  field  of  view.  Second,  a  comparison  was 

duced  by  a  small  number  of  pac  e  m  Finauy>  the  microdensitometer  tracings 

made  for  two  spot  sizes,  11  m™  “[  _  smoothed  by  eye  when  taking  the  data  from 

used  to  relate  the  spectra^  amplitudes  were  smoot  amount  of  temporal,  spatial, 

the  density  recordings.  ^gure  30  used  the  5.5-mm  spot  size  This 

and  frequency -plane  smoothing^  information,  and  the  wa\  lengths  plotted 

spot  was  selected  to  provide  higher  frequency  ^  ^  slick  is  qulte  evident  at  the 

range  from  2.9  cm  to  26  cm.  The  anomal  13.cm  plot.  Likewise,  the  over¬ 

higher  frequencies  but  is  seen  to  dm°st /^out  the  same  wavelength.  Figure  31  sub¬ 
shoot  occurring  at  1900  is  seen  to  abate  a  J^out  thejm^^  ^  ^  p.g  3Q  The 

stantiates  the  conclusions.  The  strength  than  the  spot  sizes,  the  smaller 

must  be  the  difference  in  spot  sizes.  minimUm  effects.  In  both  of  the  figures 

the  spot  size,  the  greater  the  ratio  of  maximu  limiting  the  dynamic  range  avail- 

there  is  evidence  of  the  MTF  of  the  camera-fi  y  ^  ^  diminished  as  the  analysis 
able  at  the  shorter  wa^lengths i.e  the  del ?tl i  Uh  bandwidth  of  the  interaction 

5M  -  of  data  alone. 


NRL  REPORT  7445 


31 


10  cm  fSo3 cm  Wa,elen8lh  from  .boo, 

•lized  a.  a  figure-eight  pattern  in  wavenumber  apwe  AU^fS  to  thTpomU,  **  Z 

read  along  the  major  axis  of  this  Dattprn  n..rino  a  aa  .  to  th,s  point  have  been 

~  a  abift  in  the  'oriental Jffi ttSS£E2i  SZZ?.**  ** 
measure  this  effect,  an  aperture  corresponding  to  the  wavden^mtoZ  o? lo‘to  M 

McSby^UllTSrP'X'  sfuidure  ST  T  T*  “  *■  howe, 

parent  direction  of  travel  of  the  waves  will  be  oZrwnnhT  m  Ule  ,p' 

tto\“dtor^“coutde  srrrin  the  * 

precise  **  c°M'frred  a,  indicative  of  what  happened.  Mom 

number  specttom  10  determine  the  centroid  of  the  two-dimensional  wave- 
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Figures  33  through  39  are  the  records  for  the  period  1720  to  1816  GMT  on  Oct.  5. 

The  ordinate  on  these  plots  is  not  labeled  in  dB  but  in  a  scale  of  volts,  where  IV  is  equiva¬ 
lent  x>  3  dB  of  spectral  power,  again  with  suppressed  negative  sign.  The  averaging  bandwidth 
for  'this  and  subsequent  runs  was  ±25%.  This  run  is  a  good  example  of  the  effects  of  wind 
on  the  observed  spectrum.  During  the  interval,  the  wind  was  changing  with  a  periodicity 
normally  associated  with  internal  waves.  Consequently,  the  mixing  of  variable  wind  and 
internal  waves  precludes  any  identification  of  the  individual  mechanism.  Figure  38  shows 
the  tremendous  variability  involved  in  the  8-cm  spectrum,  which  is  almost  directly  corre¬ 
lated  with  the  wind  records.  The  two  longest  wavelengths  evidence  no  such  variability. 

The  fine  structure  may  indicate  that  spectral  amplitudes  in  this  domain  are  not  greatly 
influenced  by  short-term  wind  fluctuations.  The  96-cm  curves  do  not  seem  to  change 
significantly  during  the  entire  interval,  due  in  part  to  measuring  inaccuracies  but  further 
influenced  by  the  potentially  small  number  of  wave  packets  per  unit  area  generated  by 
the  light  winds. 
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Figure  36 
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Figure  40  is  a  series  of  spectra  for  a  short  time  interval  about  1804  for  Oct.  5.  These 
data  were  obtained  by  recording  each  spectrum  on  film  and  obtaining  the  spectral  ampli¬ 
tudes  from  microdensitometer  recordings.  The  density  recording  was  averaged  in  frequency 
space  to  obtain  smoothing.  The  wind  had  been  fairly  near  zero  for  a  short  time  just  prior 
to  the  first  spectrum.  In  fact,  the  1801:45  GMT  plot  showed  no  peak  at  the  short  wave¬ 
lengths,  just  the  rise  at  the  low  frequencies.  Starting  at  the  spectra  plotted,  it  can  be  seen 
that  the  spectrum  grows  very  rapidly,  initially  at  the  highest  frequencies  and  spreading  o 
the  lower  ones  as  time  progresses.  The  curves  show  amplitudes  during  the  initial  stages 
of  growth  in  excess  of  the  dashed  line  drawn  on  each  of  the  curves.  Later,  after  equi¬ 
librium  has  been  approached,  the  spectrum  falls  and  stays  at  the  level  indicated  by  the 
dashed  line  From  these  curves  it  is  obvious  that  the  transient  spectral  response  exhibits 
and  overshoot  mechanism.  Such  an  effect  might  not  be  observable  with  a  frequency 
analysis  alone,  however,  i.e.,  a  wave  pole.  It  is  possible  that  during  the  initial  stages  of 
wave  growth,  the  spectral  energy  is  restricted  to  a  narrow  m.nge  of  angles  and  diffuses 
in  angle  with  time  until  an  equilibrium  directivity  is  achieved.  This  would  make  over¬ 
shoot  an  observable  for  directionally  sensitive  sensors  but  obscure  the  effect  for  devices 
which  measure  frequency  spectra.  Unfortunately,  the  data  collected  during  this  run  are 
insufficient  to  provide  a  definitive  answer.  The  importance  of  this  effect  on  the  phenom¬ 
enology  of  internal  wave  interactions  with  surface  waves  cannot  be  minimized.  If  the 
spectral  collapse /growth  commonly  observed  in  the  narrow  slick-like  bands  is  a  combina¬ 
tion  of  wave-wave  interaction  and  wave-growth  mechanisms,  the  interpretation  of  the 
interaction  mechanism  will  be  quite  complicated. 


WAVELENGTH  (CM) 


TIME 


Fig.  40  —  Temporal  variation  of  the  spectra  of  Oct.  5,  1971 
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Figures  41  through  47  show  the  spectral  histones  exhibited  at  the  site  during 
time  interval  2005  to  2052  on  Oct.  5.  The  film  for  this  run  was  placed  at  the  46-cm 
optical  bench  point,  indicated  as  L/2;  the  spot  size  is  5.5  mm.  The  seven  figures  showa 
very  rapid  oscillation  in  spectral  amplitude  with  a  major  dip  at  2015  and  subsequent^ 
at  about  5-min  intervals  with  diminishing  intensity.  This  is  mindful  of  a  ringing  effect 
af tersT transient  and  could  be  due  to  internal-wave  patterns  which  have  been  observed  from 
the  thermistor -chain-array  data.  As  the  water  wavelength  analyzed  approaches  ^^tes 
wavelength  the  variation  of  the  spectrum  diminishes.  This  is  probably  due  to  a  decrease 
in  the  recorded  signal  rather  than  a  real  effect.  The  small  slopes  recorded  on  film  at  these 
high  spatial  wavelengths  are  far  enough  down  on  the  modulation  transfer  curve  to  exhib 
low  signal-to-noise  ratios.  It  can  be  seen,  however,  that  the  effects  attributed  to  swell  or 
to  packet  distributions  are  still  evident,  though  of  smaller  amplitude  This  means  that 
forPfilm  recording  the  signal  levels  at  these  wavelengths  are  quite  sma  1.  In  terms  of  t 
number  of  degrees  of  freedom  contained  in  the  sampled  area,  a  data  length  of  5.5  mm 
and  a  film  spatial  frequency  of  about  20  lines/mm  for  the  highest  frequency  of  a  wave¬ 
length  of  3  cm,  the  measurement  should  be  quite  precise.  The  total  numb®"  of  d?^eeS 
of  freedom  is  about  5.5  X  20  X  20  for  the  area.  The  conclusion  here  would  be  that 
better  light  sensitivity  is  needed  to  adequately  record  such  short  wavelengths. 
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Figures  48  through  52  show  another  kind  of  effect.  These  data  were  taken  on  Oct.  6 
between  1645  and  1722  GMT.  There  is  but  one  strong  event  during  this  interval,  that  oc 
curring  about  1715.  It  is  interesting,  however,  that  the  signature  of  the  event  changes  with 
the  wavelength  analyzed.  At  the  shortest  wavelength,  i.e.,  6  cm,  the  event  shows  a  major 
dip  with  a  little  overshoot  following.  At  the  longest  wavelength  there  is  no  dip,  but  only 
a  rise  above  the  ambient  level  and  a  return  to  that  level.  Intermediate  wavelengths  show 
a  progression  from  a  dip  to  a  rise  tendency.  This  even  again  points  out  the  complexity 
of  the  interaction  being  studied.  Correlation  with  other  sensors  is  a  must  if  this  run  is  to 
provide  any  quantitative  information. 
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Figure  52 


Figures  53  through  56  represent  still  another  situation.  This  data  was  taken  on  Oct.  6 
between  1855  and  1938  GMT.  Here  the  fine  structure  is  of  large  amplitude,  and  any 
internal-wave  periodic  modulation  is  buried  in  the  noise.  Drawing  a  smooth  curve  through 
the  spectral  histories  does  yield  a  small  effect  of  perhaps  ±l-dB  modulation  at  a  period  of 
6  to  7  min,  buried  in  the  swell  fluctuation  of  ±1.5  dB.  An  obvious  interpretation  of  this 
situation  is  that  the  surface  currents  due  to  swell  were  larger  than  those  due  to  the  internal- 
wave  system.  Further,  the  wind  level  may  have  been  sufficiently  high  to  preclude  the  for¬ 
mation  of  biological  slicks  which  might  show  the  long-term  coherence  of  the  internal-wave 
system.  Waveheight  measurements  and  thermistor-chain  records  should  be  examined  for 
this  run  to  assess  the  validity  of  the  above  hypothesis. 

Figures  57  through  63  show  the  spectrum  at  seven  wavelengths  for  the  hour  1607-1707 
on  Oct.  7.  This  run  again  shows  the  frequency  dependence  of  the  observed  spectrum  and 
how  completely  different  patterns  occur  in  different  spectral  regions.  The  longest  wavelength 
shows  a  major  event  occurring  at  about  1642  in  which  only  the  level  of  the  spectrum 
changes.  At  24  cm,  one  can  see  large  changes  in  the  spectrum  in  the  interval  1612  to  1637 
GMT.  The  spectral  signal  varies  by  about  15  dB  during  that  period.  The  16-cm  graph  ex¬ 
hibits  a  behavior  indicative  of  a  pair  of  slicks,  one  at  about  1620  and  another  at  1635.  Here 
the  signal  between  the  slicks  assumes  a  level  of  about  15  dB  higher  than  in  the  slick.  Again, 
as  shorter  wavelengths  are  approached,  the  entire  level  of  the  variability  diminishes.  This 
curve  is  perhaps  the  best  evidence  obtained  for  a  frequency-dependent  interaction  of  internal 
and  wind  waves.  The  signal  at  about  1632  differs  drastically  even  in  the  range  of  8  to 
20  cm. 
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Figure  63 


Fibres  64  through  70  *.»  the  spectnim '"XK* 
This  run  is  about  the  same  as  the  Pre'=e  *I'e  attempt  to  obtain  information  at 

The  analysis  a^Krformed  m^B  ^  ^  ^  a  ,ljght  d,  „  about 

shorter  wavelengths.  In  this  case  even  a  another  diD  which  is  even  wider,  be- 

1747,  but  on  Pro“^,.n81.%‘h°rrt<!;e™eei '  ft  aiout  1740,  fte  spectrum  rises  about  8  dB 
gL"S  ■  °W  ■  -•  might  indicate  the  existence  of  two  internal-wave  crests 

Lplted”  by  about  15  min  ln 1  "f* m^hemaT^p"^^ '“ad'utUe  'or*  nolffect 

"wSngft.^  inconceivable  'that  the' 

one  may  have  had  naturally  occurring  oils  i ' X"8  a  long  period 

wave,  of  a  different  bandwidth  and  magnitude  Burn tti« '  othc”"c8j  s"“ing  the\^. 

waves  could  transfer  material  from  one  internal-wave  syste m >  t»  ^^"^e  in 

rst-ssr  sensots  wm 

have  to  be  used  to  unravel  the  questions  raised  by  this  run. 
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Figure  69 
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Figure  70 


CONCLUSIONS 

,  The  opt.^1  measurements  performed  at  the  MiiDnr*  < 

the  summer  of  1971  were  highly  successful  A  hernf?  "  San  DieB°  durinB 
directional  data  was  obtained  and  analyzed  in  the  t  ore  ,r"Posslbly  huge  amount  of 
waves  and  surface  waves.  More  than  2000  snerS  St  dy  °f  the  lntoraction  of  internal 
quency  dependence  of  this  ^  investigate  the  fre- 

features  o,  the  process  which  at  this  ^ 

that  each  ren^mTto  ™wi,U  «Me2Trffchtodire,“  7*  cxpcrlenccd-  result  is 
nlficantly  greater  body  of  data  would  need  to  ho  fcr<,|UPattrni.  It  appears  that  a  sig. 
other  parameters  to  determine  the  physte  oftoe  SSl  correla“  »*«>  sensors  of 
tion  consists  of  a  wave-generation  mnJh  f  th<?  lntcract,on  Process.  The  observed  situa- 
currenu,  surface  (Pockets),  body 

perature  effects,  and  many  other  independent !,!«  .  ’  ®  face  contami»ants,  tem- 

of  the  interaction  is  available  from  the  optical  data^Frnm  th ' Insl^ht  into  the  Physics 
be  determined  that  the  interaction  effects  are  far  if'  °pt,cal  data  alone.  it  can 

tral  fluctuations  of  up  to  15  dB  have  been  recorded  Further "th^’- nonsjationary  sPct" 
surface-wave  modulation  effects  of  ahnnt  i  rin  d6d'  Further’  there  ls  evidence  that 

perturbations  exceeding  1  dB.  It  should  be  expeLTtluhe  ***"  °Cean  SWe" 
program  could  lead  to  a  more  sophistical  mnpected  that  thc  experience  gained  in  this 

with  more  certainty  the  fundamental  principles  iSS^rSte^SST”  ‘° 
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